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ABSTRACT 

Estimation of the black hole mass in bright X-ray sources of nearby galaxies is cru- 
cial to the understanding of these systems and their formation. However, the present 
allowed black hole mass range spans five order of magnitude (1OM < M < 1O 5 M ) 
with the upper limit obtained from dynamical friction arguments. We show that the ab- 
sence of a detectable optical counterpart for some of these sources, can provide a much 
more stringent upper limit. The argument is based only on the assumption that the 
outer regions of their accretion disks is a standard one. Moreover, such optically dark 
X-ray sources cannot be foreground stars or background active galactic nuclei, and 
hence must be accreting systems residing within their host galaxies. As a demonstra- 
tion we search for candidates among the point-like X-ray sources detected with Chan- 
dra in thirteen nearby elliptical galaxies. We use a novel technique to search for faint 
optical counterparts in the HST images whereby we subtract the bright galaxy light 
based on isophotal modeling of the surface brightness. We show that for six sources 
with no detectable optical emission at the 3— sigma level, their black hole masses 
Mbh < 5OOOM . In particular, an ultra-luminous X-ray source (ULX) in NGC 4486 
has Mbh < 1244M . We discuss the potential of this method to provide stringent con- 
straints on the black hole masses, and the implications on the physical nature of these 
sources. 

Subject headings: accretion, accretion disks, galaxies:photometry, X-rays:galaxies 



1. Introduction 

Compact, off-nuclear X-ray point sources in nearby galaxies, with luminosities 10 39 — 10 41 ergs 
are referred to as Ultra- Luminous X-ray sources (ULXs). Detected in the early 1980's, with the 
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Einstein X-ray sat ellite ( Fabbianoll989 ), these objects were further studied with ROSAT ( Colbert & Mushotzky 
1999|) and ASCA (|Makishima et al.ll200d) . The XMM-Newton and Chandra X-ray observatories 
with their significa ntly higher angular resolution, dramatically confirmed the presence of ULXs 
(IKaaret et a/.ll200ll). and have enabled their spectral and temporal pr operties to be studied in detail 
(see iMiller & Colbertll2004l : iMushotzkvlEool bood: lRobertsll2007l for reviews). 



The observed luminosities of ULXs exceed the Eddington limit for a 10M Q black hole. Since 
ULX are off-nuc lear sources, their masses must be < 1O 5 M from dynamical friction arguments 
(Ka aret et al. booth . Thus, ULX may represent a class of Intermediate Mass Black Holes (IMBHs) 
whose mass range (10M© < M < 10 5 M Q ) lies b etween that of stellar m ass black holes and super- 
massive black holes observed in galaxy centers (|Makishima et al.ll2000() . Alternatively, ULX may 
be stellar mass black hole systems exhibiting s uper-Eddin gton accretions with their radiation geo- 
metrically beamed (|Shakura & Sunyaevll 19731 : Kin g 2008). 



X -ray spectroscopy has provided supporting evidence in favor of IMBHs of ~ lOOOM^ in 
ULXs dMiller et a/.ll2003l : lMiller et al.ll2003l : lcropper et al.ll2004l : bewangan et al.l2004l : [Roberts et al. 



2005; Devi et al. 



2008). Moreover, X-ray timing characteristics, i.e. presence of low frequency 

100- 



1000Mn black holes ( 


Strohmaver & Mushotzky 


2003; 


Dewangan et al. 


2006; 


Mucciarelli et al. 


2006; 


Strohmaver & Mushotzky 


2009|). While indicative these results are not conclusive, since 



there are also several arguments against IMBHs in ULXs (see e.g., iMushotzkyl 120041 : Roberts 
2007|) and further investigations are required to reveal the true nature of these sources. 



Study of the host galaxy properties of ULXs reveals that their number and total X-ray lumi- 
nosity is related to recent star formation activity , suggesting that they originate in young short-lived 
systems (lSwartz et a/.ll2004t ISwartz et al.ll2009|) . While the number of ULXs per galaxy is roughly 
the sam e for both spirals an d ellipticals, the ones in the spirals have significantly hig her lumi 



nosities (Swartz et al. 



2004|; iKuntz et al 



2005 



20041). Optical counterparts have been repor ted for some ULX (|Liu et al. 



Ramsey et all 2006 



parts have been identifie d as O stars (Liu et al 



parts are stellar clusters (|Goad et al.ll2002l : IPtak et al 



2002b 



erashima et al.ll2006). While some of the counter- 



2007), for most ULXs, the optical counter- 



2006). However, for many ULX, the opt ical 



counterparts reveal that they are either background AGN (|Gutierrezl 12006c iBonfini et al.ll2009|) or 
foreground stars. ULX s found in elliptica l galaxies may have contamination from background 
sources at ~ 44% level (|Swartz et a/.ll2004|) . Detailed stu dies of X-ray sources i n general and their 
connection with globular clusters have been undertaken kirn et alJl20o3. b009h who note that the 
X-ray properties of the the sources in the field (i.e. without optical counterparts) are not different 
from those in globular clusters. 

The allowed black hole mass range for X-ray sources in nearby galaxies span five orders of 
magnitude (1OM < M < 10 5 M Q ) and it is important to obtain tighter constrains. Here, we show 
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that the absence of a detectable optical emission allows us to impose an upper limit on the black 
hole mass for these accreting systems based on some standard assumptions. Moreoever, We argue 
that these optically dark X-ray sources cannot be foreground stars or background AGN and hence 
are a true sample of sources located within the host galaxy. To demonstrate the technique, we 
search for candiatates among X-ray bright sources detected by Chandra in archival H ST ACS, and 
WFPC2 images. 



2. Observations and Data Reduction 



Swartz et all (|2004h analyzed Chandra observations of 82 nearby galaxies and have identified 
bright X-ray s ources in them. O f these we select 13 elliptical galaxies for whi ch HST observations 
are available. Devi et al.1 (120071) analyzed a subset of thirty galaxies from the lSwartz et al\ (I2004) 
sample, and fitted the X-ray points sources with both a power-law and a disk black body model, 
hence obtaining a more conse rvative and robus t estimation of their X-ray luminosity. Nine out 
of the sub- sample analyzed by lDevi et al.1 (|2007|) are ellipticals, and for these galaxies we use the 
X-ray luminosity and coordinates giv en by them for the present analysis. For the remaining four 
galaxies we use the values quoted by ISwartz et al.\ (|2004) . The details of the sample galaxies are 
given in Table [TJ 

The optical study was carried out using the images taken with the Advanced Camera for 
Surveys (ACS), and Wide Field and Planetary Camera (WFPC2) that are available in the HST 
data archive. When observations from both cameras were available, the datasets with the longer 
exposure time, and multiple filters were preferred. For the same combination of camera and filters, 
various images available for different pointings of the same galaxy were analyzed when available, 
so as to maximize the number of X-ray sources because for different pointings, a different set of 
X-ray sources would fall in the field of view. 

Most of the optical sources are too faint to be detected against the dominant galaxy light that 
fills most of the HST images. To enhance the contrast, and aid in the detection of point sources 
in the image, the galaxy light was modeled based on the isophotes obtained using the ellipse task 
in IRAF/STSDAS software. Bright objects, if any, were masked during the fitting. The residual 
image was obtained by subtracting the model image from the observed galaxy image. The object 
extraction was done on the residual image using SEXTRACTOR with a threshold level of 3 sigma. 
On visual inspection, we find that many of the Chandra X-ray sources have counterparts in the HST 
images, within a positional offset of less than a few arc-seconds. This constant offset was applied 
for a given galaxy, to match the Chandra sources to the optical sources in the SEXTRACTOR 
catalog. When images are available in multiple filters, one of the optical images is considered as 
the reference image, and the necessary offset is applied to match it to the X-ray co-ordinates. The 
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Table 1. Sample Galaxy Properties 



Galaxy 


Distance (Mpc) 


N x 


N d 


NGC 1399 


18.3 


26 


4 


NGC 4649 


16.6 


12 


5 


NGC 4697 


11.8 


11 


3 


NGC 1291 


8.9 


5 


1 


NGC 4365 


20.9 


4 


3 


NGC 1316 


17.0 


7 


3 


NGC 4125 


24.2 


3 


3 


NGC 3379 


11.1 


3 


1 


NGC 4374 


17.4 


2 


1 


NGC 4486 


15.8 


5 


2 


NGC 4472 


15.9 


1 





NGC 1407 


17.6 


2 


2 


NGC 4552 


15.9 


3 






Note. — (1) Host galaxy name; (2) Dis- 
tance to the galaxy ; (3) Number of X- 
ray sources within HST field of view; (4) 
Number of X-ray sources without optical 
counterparts 
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(a) NGC 4486 (b) NGC 4697 (c) NGC 4649 




(d) NGC 4374 (e) NGC 1 399 (f) NGC 1316 

Fig. 1 . — : The HST images of the first six X-ray sources (marked X) listed in Table [21 Overlaid 
are I" and 3" circles centered on the shifted Chandra positions. The plus signs mark the original 
Chandra positions. Note that for four of the images, there are no optical sources even within 3" of 
the X-ray position. 
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Table 2. The Upper limit of Black hole Mass of some of the Optically Dark X-ray Sources 



Galaxy 


R.A. (72000) 


Decl. (72000) 




HST Filter 


F v x 10~ 30 


Fx/Fo 


Mu(M & ) 


NGC 4486 


12 30 50.82 


+12 25 02.66 


39 


17 +0.05 
• 1 ' -0.04 


F475W 


0.409 


533 


1244 


NGC 4697 


12 48 33.20 


-05 47 41.17 


38 


84+0-06 
• 84 -0.05 


F475W 


0.752 


243 


2890 


NGC 4649 


12 43 41.90 


+11 34 33.83 


38 


1+0-47 
yi -0.11 


F475W 


0.402 


164 


3073 


NGC 4374 


12 25 01.54 


+12 52 35.59 


39 


, +0.67 
iu -0.22 


F475W 


0.441 


347 


3378 


NGC 1399 


3 38 25.92 


-35 27 42.37 


38 


62 +0.12 
• DZ, -0.09 


F606W 


0.370 


228 


3927 


NGC 1316 


3 22 36.46 


-37 13 24.68 


38 


on+0.23 
• 8U -0.08 


F475W 


0.449 


179 


4780 


NGC 1316 


3 22 35.58 


-37 13 14.10 


38 


78+0.40 
' 70 -0.07 


F555W 


0.520 


287 


6366 


NGC 1399 


3 38 32.33 


-35 26 45.73 


38 


C4+1.13 
• J ^-0.41 


F606W 


0.377 


186 


7829 


NGC 1399 


3 38 27.62 


-35 26 48.76 


39 


49+0.16 
•^ z -0.14 


F606W 


0.766 


702 


7829 


NGC 4649 


12 43 34.17 


+11 33 41.93 


39 


04+0.10 
u ^-0.11 


F475W 


0.912 


97 


8073 



Note. — (1) Host galaxy name; (2) Right Ascension of shifted position in hours, minutes and 
seconds; (3) Declination of shifted position in degrees, arcminutes and arcseconds; (4) log of X-ray 
luminosity in ergs/s; (5) HST filter for which the upper limit on flux and the black hole mass limit 
is calculated; (6) Upper limit on Optical flux in ergs/s/cm 2 /Hz; (7) Lower limit on ratio of X-ray to 
optical flux; (8) Mjj, Upper limit on black hole mass. 
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images in other optical filters are then aligned to the reference image using geomap and geotran 
tasks in IRAK While a more detailed report on the nature of sources with optical counterparts will 
be presented later, in this work, we concentrate on those sources for which no optical counterpart 
was detected. 

X-ray sources which did not have an optical counterpart (at the 3— sigma level) within I" of 
their shifted positions, are termed as optically "dark" sources. Six examples of such sources are 
shown in Figured! where the I" and 3" circles centered on the X-ray co-ordinates are overlaid on 
the observed optical image. Note that for four of these images, there are no optical source within 
3" of the X-ray position. Having identified such optically "dark" X-ray sources, we estimate the 
upper limit on their optical flux based on the 3— sigma threshold at that position. 



3. Optically dark X-ray sources 

These optically dark sources are X-ray bright compared to their optical emission and hence 
are not foreground stars. This can be further quantified by estimating the X-ray-to-optical flux 
ratio log (fx / fo) where fx is the unabsorbed flux in the 0.3 — 8 keV band and fo is the flux in 
an optical band. This ratio ranges from 0.1 to 50 for AG Ns including BL L acs when fx is in 
the 0.3 — 3.5 keV band and the V-band magnitude is used (|stocke et al. 1991 ). In contrast, the 



estimated lower limit for the optically dark sources in the sample is significantly larger. This is 
illustrated in Table © where the ratio is given for ten sources. We, therefore, conclude that these 
sources are not background AGNs. 

Thus, these sources are most likely to be bright X-ray binaries (or at least accreting systems) 
within the galaxy. An accretion disk around a compact object should also prod uce optical emission 



whos e flux can be estimated as follows. In the standard accretion disk theory (Sh akura & Sunyaev 



19731) . the effective temperature profile as a function of radius R, of an accretion disk around a 



black hole with mass M and accretion rate M, is given by oT 4 (R) = ^ G! ^i M 8(R) where 8(R) — 
1 — (6GM/Rc 2 ) 1//2 . The observed flux from the disk at a frequency v is then given by the integrated 
sum of the black body emission over all radii, 

cos? r R ° u ' , , NN 
p v = -RT B v (v,T(R))2nRdR (1) 

where B v (v,T(R)) is the blackbody intensity, R out and = 6GM/c 2 are outer and inner radii of 
the disk, i is the inclination angle of the disk and D is the distance to the source. Assuming that 
most of the contribution to F v arises from regions in the disk that are far away from the inner and 
outer radii, the expected observed flux can be written as 

F V ~ 7xl0- 31 ergss- 1 cm- 2 Hz- 1 (^)- 1 / 3 (i? r )- 2 / 3 
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X f h. D ) -2 ( M N2/3 (2) 

l 10 39 ergS5- lJ l 10Mpc J 1 lOOOMg/ V } 

where A is the wavelength and r\ = L x /Mc 2 is the radiative efficiency of the accreting system and 
cos/ is taken to be 0.5. 

For optically dark ULX, the predicted accretion flux should be less than the measured upper 
limit F V)max . Thus one can estimate an upper limit on the black hole mass as 

Mu < 1000 M ( Fv >- 1 ) 3 / 2 

7 x I0~ il ergs s _1 cm~ 2 Hz 1 

y( ^ \VV 11 )( Lx D )3 (3) 

V 5000A ; v 0.1 n 10 39 ergs^ i; v 10Mpc ; v ' 

For each dark X-ray source in our sample, and for all available filters, we estimate this upper limit 
on the black hole mass. We use the integration (Eqn Q]) to evaluate the upper limit, rather than 
the approximation Eqn. §3$ i.e. we take into account the effect of the inner boundary condition 
(8(R) = 1 — (6GM/Rc 2 ) 1//2 ) on the temperature profile. The difference in the upper limit obtained 
is marginal (< 20%). We assume a standard radiative efficiency of r\ = 0.1 and cosz = 0.5. To 
obtain a more conservative upper limit, the one-sigma lower value of the X-ray luminosity are 
used L x — AL X . In Tabled we list the ten best cases in ascending order of black hole mass limit 
Mjj. For the other sources, My > 10000M Q and hence is not a significant constraint. The best case 
is for the ULX in NGC 4486, for which M v = 1244M©. 

X-ray irradiation of the outer disk may increase the local temperature there and the disk may 
emit a large r optical emi s sion. We have estimated this effect using the formalism given in the 



appendix of IVrtilek et al.l (119901) and find that X-ray irradiation is not important for the constraint 
obtained here. 



4. Discussion 

Optically dark X-ray sources cannot be foreground stars or background AGN, otherwise their 
optical emission would be significantly higher than what is detected. Hence, these are a clean 
sample of sources within the host galaxies, which are probably accreting black hole systems. The 
optical emission from a standard accretion disk scales as mass of the black hole M 2 / 3 and hence 
the non-detection of optical emission imposes an upper limit on the black holes mass Mu- For 
ten of the sources Mjj < 10000M©. For a source in NGC 4486 with an X-ray luminosity clearly 
exceeding 10 39 ergs/s (and therefore a bona-fide ULX by definition), the estimated black hole mass 
is smaller than 1244M Q . This is two orders of magnitude smaller than the constraint obtained from 
dynamical friction, which is 1O 5 M . 
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These sources with black hole mass, Mjj < 5000M Q cannot be accreting systems with massive 
black holes residing in star clusters, or i n the nuclei of merged satellite galaxies. For typical low- 



luminosity dwarf galaxies (Mb ~ — 8.0: lMatedll998[) . such an optical counterpart would be easily 



detected, given that our 3— sigma limits on the HST images are much fainter. Even a compact 



nucleus of a merged dwarf galaxy (Mg ~ —7.0; lLotz et al.ll2004|) would have been easily identi 



fled. If they are binary systems, their companion cannot be a massive O star as such a star would 
have been detected in the optical image. Assuming an O star, with Mg ~ —5.5, we find that the 
possibility of such a companion can be ruled out in all cases for which M\j < 5OOOM (Table©. 

In all of the above arguments, we have ignored the effect of dust obscuration in the host 
galaxies, because we are only considering elliptical galaxies in the present work. Although some 
ellipticals are known to have dust lanes, and nuclear rings in their centers, the X-ray sources we 
are considering here are distributed at fairly large radial distance from the center to be significantly 
affected by dust. 

Even for the best cases of optically dark X-ray sources presented in Table ©, the range of 
black hole mass allowed is still large as the source could be a & 1000A/O intermediate mass black 
hole, or a few solar mass object emitting at super-Eddington luminosities. The brightest X-ray 
sources in the sample have a luminosity of a few times 10 39 ergs/s. This is unfortunate, since a 
dark source with luminosity > 10 40 would have provided an order of magnitude better constraint 
on the black hole mass. Since the black hole mass upper limit My «D 3 a bright X-ray source 
in a more nearby galaxy (D ~ 2Mpc) would have also provided significantly better constraints. 
A systematic search for such sources in very nearby galaxies may indeed prove fruitful. Another 
point to note is that bright X-ray sources in these galaxies are known to be variable in X-rays. 
Our analysis in this work, implicitly assumes that the X-ray luminosity observed through a single 
Chandra observation, represents an average luminosity which is used to derive an average accretion 
rate M = Lx/ (tic 2 ) which in turn is used to estimate the upper limit on the black hole mass (Eqn 
|3]). This assumption is required because the expected optical emission arises from the outer part of 
the disk and the local accretion rate there may be different than the one in the inner region which 
produces the X-rays. Any accretion rate fluctuation in the outer disk will be transfered along the 
disk on the viscous time-scale which could be significantly longer than a day. Thus, in principle 
one needs to ascertain the average X-ray luminosity of a source, and using a single very bright, but 
rare, X-ray observation of the source will not represent the average accretion rate. 

A systematic and comprehensive multi-wavelength study (using also other bands like infra- 
red and radio), along with X-ray variability studies, can shed further light on the nature of these 
sources. 
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